Perturbation theory for plasmonic eigenvalues 



Daniel Grieser and Hannes Uecker 
Institut fur Mathematik, Carl von Ossietzky Universitat, D-26111 Oldenburg, Germany 



o 



o 
o 



> 

(N 

in 

rn 
O 



X 



Svend-Age Biehs*, Oliver Huth, Felix Riiting, and Martin Holthaus 
Institut fur Physik, Carl von Ossietzky Universitat, D-26111 Oldenburg, Germany 

(Dated: October 27, 2009) 

We develop a perturbative approach for calculating, within the quasistatic approximation, the 
shift of surface resonances in response to a deformation of a dielectric volume. Our strategy is 
based on the conversion of the homogeneous system for the potential which determines the plas- 
monic eigenvalues into an inhomogeneous system for the potential's derivative with respect to the 
deformation strength, and on the exploitation of the corresponding compatibility condition. The 
resulting general expression for the first-order shift is verified for two explicitly solvable cases, and 
for a realistic example of a deformed nanosphere. It can be used for scanning the huge parameter 
space of possible shape fluctuations with only quite small computational effort. 
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I. INTRODUCTION 



Recent advances in plasmonics, i.e., in the use of sur- 
face plasmons for subwavelength optics^— have led to 
renewed interest in the physics of plasmon excitations 
bound to smooth or rough surfaces^ In particular, it 
has been proposed to transport electromagnetic energy 
along linear chains of nanoparticles, 15 possibly embedded 
in a gain medium^ and proof-of-principle experiments 
have been madei^r— In such setups, fabrication-induced 
shape imperfections of the nanoparticles inevitably will 
result in slight shifts of their resonance frequencies, and 
it might be crucial to estimate the maximum size of such 
fluctuations which can still be tolerated*^ Although so- 
phisticated numerical methods for computing the optical 
response of nanoparticles do exist ) 11 ' 12 it would still be 
helpful to have a flexible analytical tool which exploits 
the fact that the shape of an unintentionally deformed 
nanoparticle is close to some theoretical ideal, as this tool 
would allow one to explore the huge parameter space of 
possible perturbations in a computationally cheap man- 
ner. 

In the present paper we develop such a perturbative 
approach to the computation of surface resonances. Our 
strategy is quite general, relying on concepts borrowed 
from differential geometry^ The mathematical argu- 
ments are given in the following Sec. [TTJ some technical 
details have been deferred to Appendix [AJ The main 
results are the expressions (f26|) and (l28l) . which quan- 
tify the shifts of the resonance frequencies to first order 
in the deformation strength. Section IIIII then provides 
three illustrative examples. The first two of these make 



contact with analytically solvable models, thus helping to 
gain confidence in the formalism, while the third one is a 
more realistic application to a deformed nanosphere for 
which no closed analytical solution is available. This ne- 
cessitates to consider the splitting of degenerate modes, 
which is done in Appendix [Bj The paper ends with some 
concluding remarks in Sec. IIVI 



II. 



PERTURBATION THEORY FOR 
PLASMONICS 



Consider a volume O C I 3 filled with a dielec- 
tric medium which is characterized by an isotropic, 
frequency-dependent dielectric function e(w); this volume 
fl be surrounded by vacuum. We employ the quasistatic 
approximation, which is valid if the characteristic linear 
extensions of are small in comparison with the wave- 
lengths of impinging radiation, and thus captures essen- 
tial features of nonretarded plasmon dynamics in small 
nanoparticleSi 12 ' 14 " — The Fourier components u(r\u>) of 
the potential then are given by the solutions to the set 
of equations 



Am = in K \ dil 
w_ — u + = on dQ 
ed n U- — d n u + = on dQ . 



(1) 
(2) 
(3) 
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Here U- is the restriction of u to the interior VL and u + its 
restriction to the open exterior M 3 \ f2; both it_ and u+ 
are smooth everywhere except at the boundary d£l. We 
stipulate that dQ. be sufficiently smooth to substantiate 
the following operations, and suppress the dependence of 
u on the frequency ui altogether. The expressions d n u± = 
n • Vw± in Eq. ^ are the derivatives of the potential in 
the direction of the outward unit normal n. 

Observe that at this point the knowledge of the dielec- 
tric function e(u>) is not yet required. Rather, we take e 
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on the l.h.s. of Eq. ([3]) as a real number and regard it as an 
eigenvalue, henceforth dubbed as plasmonic eigenvalue, 
if the system (|TJ) — (J31) with that particular e possesses 
a nontrivial solution u which vanishes at least quadrati- 
cally with increasing distance from dVL. This nomencla- 
ture is formally justified by reformulating the system ([1]) 
- ([3]) in terms of Dirichlet-to-Neumann operators, so that 
the desired values of e explicitly appear as inverse eigen- 
values of a combination of such operators. 18 This some- 
what unfamiliar view offers conceptual advantages: On 
the one hand, the material-specific aspects of the prob- 
lem, embodied in the function e(uj), are separated from 
the geometric ones; only these geometric aspects matter 
when studying deformations of fl. On the other hand, 
standard theorems concerning the behavior of eigenval- 
ues of linear operators under perturbations can now be 
employed. After such a plasmonic eigenvalue has been 
found for a particular geometry, i.e., for some domain Q 
with a given shape, that eigenvalue then determines the 
frequency w s of a surface resonance through the equation 
e = Re(e(w s )). For example, a half space with planar 
boundary dfl yields e = — 1, and the dielectric function 
can be taken as 
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(4) 



for metallic materials with plasma frequency w p , describ- 
ing the free electron motion. This then leads to the fa- 
miliar expression uj s = cj p /v2 for the surface plasmon 
resonance at a planar metal- vacuum interface in the qua- 
sistatic limits Typical plasma frequencies for good con- 
ductors are on the order of 10 16 s _1 . 

We now consider the response of a plasmonic eigen- 
value e to a small deformation of ft. This deformation is 
modeled in terms of a shape function a : dQ — > M and a 
dimensionless parameter h, such that the surface of the 
deformed volume is given by 



dft(h) = {r + ha(r) n(r) : r <= dfl} . 



(5) 



We assume that one has solutions e(h), u(r,h) to the 
problem ([J) - ([3]) with fl replaced by Sl(h), and that 
e(h) and u(r, h) depend diffcrcntiably on h. We denote 
the derivatives in h by a dot and write e = e(0) and 
u(r) — u(r,0) for brevity. 

In what follows we derive a system of equations for e 
and u(r), and therefrom an explicit expression for e. To 
this end, we take the derivative of the system (JTJ - © 
with respect to h. First, for any r g" dfl = dfl(0) one 
has r dfl(h) for h sufficiently close to zero; therefore, 
Eq. ([!]) readily yields 



Aii = in R 3 \ dn . 



(6) 



Next we differentiate Eq. @. In a more explicit form, 
this equation reads 

it_(r + ha(r) n(r), h) — u + (r + ha(r) n(r), h) = (7) 



and hence gives 

an • Vu_ + u_ — (an • Vw,+ + u + ) = (8) 



u + = an • Vu+ — an • Vm_ 
= (e - l)ad n u- , 



(9) 



where Eq. ([3]) has been used for eliminating u + . 

Similarly we differentiate Eq. ([3]) , keeping in mind that 
d n u± — n • Vti±, and taking into account that n also 
depends on h. Thus, one has 



(n • Vu±J — n • Vii± + n • (Vu±) 



(10) 



The required derivatives n and n • (Vu±)' are calculated 
in Appendix El According to Eq. (|A12|) . n = —Vga co- 
incides with the negative gradient of the shape function 
inside the surface dfl, while Eq. (|A21[) expresses n- (Vu±)' 
in terms of the Laplace-Beltrami operator = div^ 
applied to u±, and the mean curvature H of <9f2. Com- 
bining these results, we find 

(n • Vu±)' 

= — V^a • VgM± — aAgu± + 2aHd n u± + d n u± 
= — divg(aVgu±) + 2aHd n u± + d n u± . (11) 

Thus, taking the /i-derivative of Eq. © leads to 

(e<9„u_ - d n u + ) = ed n u- + e{d n U-) — {d n u + ) 
= ed n ii- 

— diva (a(eVau_ — Vgu + )) 
+2aH(ed n U- — d n u+) 
+ed n ii- - d n u+ . (12) 

The l.h.s. now vanishes identically, and the third term 
on the r.h.s. vanishes because of Eq. ([3]) itself. Moreover, 
since u+ — u_ on dVl according to Eq. ^j, and Vg differ- 
entiates inside dQ only, we can replace Vgu+ by Vgu_ 
in the second term, finally leaving us with 

ed n ii- — d n ii + = —ed n u_ + (e — 1) divg(aVgu_) . (13) 

In summary, the differentiation of the homogeneous sys- 
tem (P) - ([3J has led to the inhomogeneous system 

Aii = inR 3 \dn (14) 
u_ — ii+ = F on dfl (15) 
ed n U- — d n ii+ = —ed n U- + G on <9£1 , (16) 



where 



F = (e — l)ad n U- 

G = (e — 1) diva(aVa<tt_) 



(17) 
(18) 



Since the homogeneous system has a nontrivial solution 
(the given u) , its inhomogeneous descendant can admit a 
solution only if the right hand sides F, —ed n ii-+G satisfy 
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a certain compatibility condition. Since by assumption 
there is a solution (the /i-derivative of the given family 
u(r, h)), this condition actually is fulfilled, and yields an 
expression for e. For deriving this expression, we start 
from Au- — Au_ = in fi. Green's formula then gives 

= / d 3 r (u_Aii_ - (Au_)ii_) 
Jn 

= f d 2 S (u^d n ii^ - {d n u.)iiJ) , (19) 
Jan 

where d 2 S denotes the surface area element. Thus, one 
has 



d 2 Su-d n u-= d 2 S(d n u-)u- ; (20) 
an Jan 



a corresponding identity holds for u + , u + . Therefore, 
d 2 S (u_e9„ii_ — u + d n u + ) 



on 



on 



d 2 S (e(d n uJ)ii- - (d n u+)u+) . (21) 



Since u_ = u+ on 9f2, the l.h.s. becomes 

/ d 2 ^ u- (ednU- — d n u+) = / d 2 S u~ (— ed n u-+G) . 
Jan Jan 

(22) 

On the other hand, since d n u+ = e<9„u_ on 9f2, the r.h.s. 
of Eq. (|2T1) takes the form 



d^S* e(d„u_) («_ - u+) = e / d z S (<9„u_)F . 
an Jan 

(23) 

Equating these two expressions (|22|) and ([23]) , and solving 
for e, we obtain 

■ J m d 2 Su-G-ef dn d 2 S(d n u„)F 
I d n d2Su - d ^ u - 

Inserting the formulas (fT7| and ([T8|) for F and G, and 
performing an integration by parts, the numerator re- 
duces to 



/ d 2 ^ u_(e — 1) diva(aVgu_) 
Jan 



-e / d 2 S {d n u_)(e- l)ad n u_ 
Jan 

= (1-e) / d 2 Sa((V 9U _) 2 + e(cU_) 2 ) . (25) 
Jan 

Thus, the desired expression for the first-order change 
of the plasmonic eigenvalue e with the deformation 
strength h finally reads 



e = i - 



Jan d s u-dnu- 



This is the principal result of the present work. Observe 
that Au_ = implies 

/ d 2 Su_<9„7i_ = / d 3 r(V«_) 2 , (27) 
Jan Jn 

so that the denominator is positive. Of course, e can 
likewise be expressed entirely in terms of u+: 



e = U - 



f dn d 2 Sa(e(X7 d u + ) 2 + (d n u + ) 2 ) 
Jan d 2 Su+d n u+ 



III. APPLICATIONS 

Consider an infinite half-space geometry with the di- 
electric medium filling the volume z < 0, so that its 
boundary is given by the plane z = 0. Let k = (k x ,k y ) 

be a two-dimensional wave vector, with k = \Jk% + fc 2 . 

Solutions to Laplace's equation (TTJ) which vanish for 
z — > ±oo then are given by 



u+(r) 



d 2 K 



(27T) 2 
d 2 K 

(2< 



A(K)e IKX " KZ (29) 
2 B(k) c iK - x +^ , (30) 



with x = (x, y) and r = (x, y, z), and the continuity con- 
dition ^ yields B(n) = A(k). Moreover, since u± is 
real, one has A*(n) = A(—k). At this point, the ampli- 
tudes A(k) describing the excitations at an exactly pla- 
nar surface are arbitrary. The normal and the in-plane 
derivatives at z — then are 



d n u± = T 



V a u± 



d 2 K 



(2tt) 2 

d 2 K 



A(k) k e if 



(31) 



(2tt)3 



A(k) (in) e IK X , (32) 



respectively, so that the condition Q immediately pro- 
vides the well known eigenvalue e = — 1 for this particular 
geometry. We write 



(V aU ±) 2 



together with 



(2tt) 2 7 (2tt) 2 
d 2 «; /■ dV 



(2tt) 2 J (2tt) 2 



fWv(«-«')-x 
A(k)A*(k') (k-k') e i(K - K '>* 



(2tt) 2 y (2tt) 



giving 



d 2 K 



d 2 Su±cU± = T / 7^ I^WI 2 ^ ■ (33) 
an J (2tt) 2 
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Now we introduce a small deformation of the planar sur- 
face, described by some suitable shape function a(x) 
The required amplitudes A(k) consequently are deter- 
mined by that deformation; the assumption that the ex- 
act potential u±(r) can still be written in the form ([29] 
or ([30]) constitutes the Rayleigh hypothesis^ Upon in- 
serting the Fourier transform of the shape function, i.e. 



o(x) = 



d 2 K 



o(k) e ih 



(34) 



und using e = — 1 for the unperturbed eigenvalue, either 
Eq. (|26p or its variant (|28p readily yields 

. _ 2 f_wj* °(« ; - *)A(«)A»(k>«'(k • ft' - 1) 

(35) 

where & denotes the unit vector in the direction of k. 
In particular, if a(x) = c, one has a(«) = (2it) 2 c5(k), 
so that e = 0: The plasmonic eigenvalue does not 
change when the entire surface plane is displaced by some 
amount c. 

On the other hand, there exists an exact integral 
equation for determining the amplitudes A(k) associated 
with a given surface deformation, derived by Farias and 
Maradudin on the basis of the Rayleigh hypothesis^ 



e+1 
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+ 1 f d K 

A(k') = / , n .„ J(k - k\k' - k)(1 - ft ■ k')kA(k) 

- 1 J (2ny 

(36) 



with 



J(a|q) = i / d 2 S , e^ q x (exp(a/ia(x)) - l) ; (37) 

this has been applied by Maradudin and Visscher to the 
study of particular perturbations of planar surfaces. 20 
Expanding the latter expression (f3"T|) to first order in h 
gives 



J{a\ci) = ha{c L )+0{h 2 ) , 
so that Eq. (|36|) becomes 
e + 1 



(38) 



1 fan 

-A(K') = h J -j—y a(K' -k)(1-k-k')kA(k) (39) 

for sufficiently weak perturbations. Recalling that the 
unperturbed eigenvalue is —1, we insert e = — 1 + 0(h) 
into the denominator on the l.h.s. Multiplying both sides 
by A*(k')k' , integrating, and rearranging then yields 



(40) 



e = -I + he + 0(h 2 ) , 



with e indeed formally equal to the previous Eq. Q35p. 
Thus, our perturbative result is consistent with the for- 
mula ([36]). 

A second example which allows one to confirm the va- 
lidity of the perturbative approach by analytical means 
is provided by the deformation of a dielectric sphere of 



radius R into a spheroid. In this case the unperturbed 
potential is written as 

oo 

u + (r) = YA t r-^Y^) 



«_(r) = VVrJw 



£=0 



(41) 



valid for r > R and r < R, respectively. Here denote 
the familiar spherical harmonics. The restriction to these 
basis functions with m — 0, which do not depend on 
the azimuthal angle cp, confines us to deformations which 
preserve the rotation symmetry around the z-axis. From 
condition © one gets B t = A e R~( 2l!+1 \ and Eq. © then 
leads to the eigenvalues 



+ 1 



(42) 



In particular, e\ = —2 characterizes the dipole reso- 
nance^ 

We deform the sphere into a spheroid oriented along 
the z-axis, employing the shape function 



a($,ip) = i?cos 2 (i?) 



(43) 



Here the choice of the sphere's radius R as the scale of 
the deformation is a matter of convenience. The three 
degenerate dipole modes of the sphere are shifted when 
the deformation strength h adopts nonzero values, such 
that the mode associated with the z-axis splits off from 
the two others. For estimating the corresponding change 
of ei for small h, we use 



ui_(r) = rF 1 (?9) = r\ — cos$ , 

V 47T 



(44) 



and calculate 



3 

(<9„ui_) 2 = — — cos 2 1? 

47T 

(V 9Wl _) 2 = A sin 2tf 
4"7T 

3 

Ui_9 n ui_ = — Rcos 2 ■& 

4-7T 



on d£l, giving 



crSaidnm-Y = -R A 

/ d 2 Sa(V dUl _) 2 = ^R 3 
Jon 5 



d^Ui-Snlii- 



R 3 . 



on 



Plugged into Eq. (j26|) . this yields 



ei 
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(45) 



(46) 
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This easily obtainable result can be checked against the 
exactly known expression for the dipolar z-mode of a 
spheroid: 21 



1 

£1 = 1 -I 



(47) 



where L denotes the so-called depolarization coefficient, 
given by^ 2 . 



L(e) 



i^f- (artanh(e) - e) , h > 
Mf^(e - arctan(e)) , h< 



with 



(48) 



(49) 



Thus, e 2 becomes small with h, allowing us to expand 
L(e) as 



He) 



±-^e 2 + 0(e 4 ), /i>0 



3 ' I 5 



i + ^e 2 + G(e 4 ), /i<0 



(50) 



Exploiting Eq. (|49|) one expresses L in terms of h, finding 

(51) 



This tells us that ei, as given by Eq. (|4"7|) . should be 
expanded in powers of (X — 1/3): 



(52) 




Inserting Eq. ([51)1. we obtain 



e a = -2-^/i + 0(/i 2 ) , 
5 



(53) 



in accordance with the result ([46)1 of perturbation theory. 

The following third example concerns a deformed 
nanosphere which does not admit an analytical solution 
in closed form, so that the accuracy achieved by first- 
order perturbation theory has to be ascertained by com- 
parison with numerical calculations. This is exactly the 
type of application we have in mind, since here Eqs. ()26[) 
and ([2"8]) provide a quick and reliable estimate of the pos- 
sibly detrimental consequences of geometrical imperfec- 
tions^ without hard requirements on computational re- 
sources. We assume that the ideal sphere is distorted by 
two Gaussian protrusions, and parametrize its surface as 



i(i9, f ) = ft exp 



(54) 



with dist(i?j, <fi] ®, ip) denoting the dimensionless Eu- 
clidean distance between the two points on the unit 



N 







-1 




-1 
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1.1 



1.0 







-1 







FIG. 1: (Color online) Sketch of the particular perturbed 
nanosphere (|54p employed as example for demonstrating the 
accuracy of the perturbative approach. The deviation of the 
scaled surface radii r{&, <p)/R from their ideal value 1 is visu- 
alized here for h = 0.2. 



sphere specified by the angles fiijtpi and ip, respec- 
tively. The protrusions' parameters are chosen, some- 
what arbitrarily, as fi = 1.0, f% — 0.8 for their ampli- 
tudes, w\ = 0.7, u>2 = 0.65 for their Gaussian widths, 
and i?i = 1.1, d 2 = 2.2 and cpi = 4.9, <p 2 = 0.1 for their 
locations. Figured] depicts the resulting deformed sphere 
for the value h = 0.2 of the overall perturbation strength. 



This example also illustrates a further important issue. 
While the dipole modes of a perfect sphere are threefold 
degenerate, this degeneracy is removed entirely by the 
distortion ([541 . Thus, there now are three branches of 
dipole-like eigenvalues, and our previous analysis applies 
to each branch separately, provided the starting point is 
chosen appropriately: The expressions ([2f))) and (|28p re- 
fer to the individual branches, if the proper linear com- 
binations of the unperturbed degenerate modes are in- 
serted. The problem how to find these linear combina- 
tions is solved in Appendix [Bj Basically, the numera- 
tor of the formula ([2"rJ[) defines a quadratic form of the 
eigenmodes (see Eq. (jBll[) h the required proper linear 
combinations of the unperturbed modes are those which 
diagonalize this form. Their determination, and the eval- 
uation of the ensuing surface integrals determining e for 
all three branches, does not demand much numerical ef- 
fort. Figure O shows the perturbative results for defor- 
mation strengths —0.25 < h < +0.25 in comparison with 
data obtained by nonperturbative numerical computa- 
tions. 2 '' Quite remarkably, first-order perturbation theory 
still produces excellent results when the shape variation 
already is quite substantial, that is, for values of h up to 
0.1; even for h as large as 0.2 one obtains good estimates. 
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-0.2 0.2 

h 

FIG. 2: Numerically computed plasmonic eigenvalues for the 
dipole modes of the deformed nanosphere depicted in Fig. [T] 
(crosses), in comparison with the results of first-order pertur- 
bation theory (full lines). 



IV. CONCLUDING REMARKS 

The first-order perturbative expression (|26l) or (12511 for 
the shift of a plasmonic eigenvalue clearly has a limited 
range of applicability, insofar as it is restricted to suffi- 
ciently weak deformations, but it stands out because of 
its generality, and easy use. Together with similar results 
for higher derivatives^ one obtains a formal perturba- 
tion series 

e(h) ^ e + he+^h 2 e + ... . (55) 

Above we have simply assumed the existence of e(h) and 
u(r, h). Given a solution at ft. = 0, this existence actually 
follows from the reformulation of the problem in terms 
of Dirichlet-to-Neumann operators,— and from standard 
perturbation theory for eigenvalues. If e is p-fold degener- 
ate (as in the case of the sphere), then one has p branches 
€Y>, vy'\ i — 1, . . . ,p. The analysis of Sec. [TT] then ap- 
plies to each branch separately, with being given by 
Eq. (HU) or (gFJ with u = u®. Here uW form a 

basis of the subspace of degenerate solutions at h = 0, de- 
termined such that they diagonalize the quadratic form 
given by Eq. (|B11|) . 

The three examples we have given in Sec. IIIII vary in 
character, the first two recovering known analytical re- 
sults and thus confirming the correctness of our formal 
line of reasoning. The third example, summarized by 
Figs. [1] and [2J demonstrates the utility of our approach for 
practical purposes. Here we have dealt with an asymmet- 
rically deformed nanosphere, and obtained fairly good es- 
timates for the shifted dipole modes. If we assume that 
tolerances on the order of 5% can be met in nanosphere 
fabrication, first-order perturbation theory thus allows 
one to quantify the effects of a large variety of possible 
shape fluctuations with both sufficent accuracy and only 
small numerical effort. 
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Appendix A: Auxiliary calculations 

In this appendix we compute the derivatives entering 
into the evaluation of Eq. (fTU)) en route to the important 
result (|13l) . This calculation invokes some notions from 
differential geometry— We fix a point ro G <9f2 = 951(0) 
and parametrize d£l near ro as r = r(v,w), with ro = 
r(0,0). It is assumed that the parametrization is such 
that the tangent vectors r v = d v r and r w — d w r are 
orthogonal, and have unit length at ro, and such that the 
v-line (i.e., v r(v, 0)) and the w-Mae through ro have 
normal vectors parallel to n at ro (i.e., have vanishing 
geodesic curvature). 

The curvatures of the u-line and the ui-line are Ky — 
r vv ■ n and k w — r ww ■ n, respectively. Considering n as 
a function of v, w, one has 

r v ■ n = (Al) 

and hence 

r vv • n + r„ • n„ = , (A2) 

giving 

Kv — r^ • ; (A3) 

analogously n w = —r w ■ n w . 

Note that orthonormality of r„ and r w holds only at 
ro, not at nearby points. In the following calculations we 
always evaluate at h — after differentiating in h, and 
at r = ro after differentiating in v, w. 

We now derive an expression for n. The surface d£l(h) 
is parametrized by 

r(v, w, h) = r(v, w) + ha(v, w)n(v, w) , (A4) 

where we abuse notation by considering a and n as func- 
tions of v, w. Its unit normal is n = N/||N||, where 
N = r„ x r tt „ At r one has n = N by assumption. 
Differentiating, one gets 



||N|| ||N||Hl|N|| J 
= -j- (N - n(n • N)) , (A5) 

which states that h is determined by the tangential part 
of N. Moreover, 

N = r„ x r w + r v x r w . (A6) 
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Since r = an, one has 



a„n + an,, 



(A7) 



At 1*0, the identity n = r v x r w implies n x r w — —r v and 
r v x n = — r w . In addition, n v x r w — (n„ ■ r v )r v x r w , 
because the component of n„ parallel to r w drops out of 
the vectorial product. (Observe that n • n = 1 implies 
n v ■ n = 0, so that n„ lies in the span of r v and r w .) 
Putting this together and using Eq. (|A3|) . one finds 



r„ x r 1t 



(a v ii + an.„) x 

L W 

-a v r v + a(n v ■ r„)n 



analogously, 



(A8) 



(A9) 



Inserting Eqs. (|A8[) and (IA9[) into Eq. (IA6|) . we arrive at 

N = - (a v r v + a w r w ) - a (k v + n w ) n . (A10) 

Now we observe that 

a v r v + a w r w = Vga (All) 

is the gradient of a as a function in the surface dQ. Since 
according to Eq. (|A5[) the desired derivative n is the tan- 
gential component of N, we have the compact result 



h = — V«a 



(A12) 



Next we calculate n • (Vu)' for u = u± or u = u_ at ro. 
Starting from 



(Vu)'(r) 



d_ 

dh 



(Vu)(r + ha(r)n(r), h) 
a(n-V)Vtt(r) + Vii(r) 



(A13) 



one finds 



n • (Vu)' = a(n ■ V) 2 u + d n u 



(A14) 



To simplify notation we assume that coordinates in ffi 3 
are chosen such that n = (0, 0, 1) at ro, reducing the first 
term on the r.h.s. to ad^. 3 u at ro. Because Au = 0, we 
then have, at ro, 



n • (Vu)' = -a (d1 x + dl 2 ) u + d n u 



(A15) 



This has to be expressed in terms of derivatives of u inside 
the surface. We choose the surface coordinates as v = X\, 
w = X2, but still write r = r(v,w). Differentiating u 
twice with respect to v, 

d v u(r(v,w)) = (r v ■ V)u(r(v, w)) 
dvu(r(v,w)) — (r v ■ V) 2 u(r(v, w)) 

+ (r vv ■ V)u(r(v,w)) , (A16) 



and observing r v ■ V = d Xl and r^^ = K v n, we get 

d 2 v u = d 2 Xi u + Kv d n u . (A17) 
Combining this with the analogous equation for X2 gives 

(d 2 xi + O u = ( d v + d D u - ( K - + K -) 9 « u • ( A18 ) 

Introducing the mean curvature H of dQ, 



and identifying 



H — — (k v + k w ) 



A d U = (dl + d 2 w ) U 



(A19) 



(A20) 



as the Laplace-Beltrami operator of the surface dtt, 
Eq. (|A15|) takes the final form 



n • (Vu)' = — aAgu + 2aHd n u + d n u . 



(A21) 



While we had used a special coordinate system for sim- 
plifying the derivation, the Laplace-Beltrami operator 
is defined independent of the choice of coordinates as 
Ag = div gVg. Here div^ is the divergence inside the 
surface, given by the negative adjoint of the gradient with 
respect to the surface volume element. This is exploited 
in Eq. {nj. 

The fact that the two results (|A12|) and (|A21|) obtained 
in this appendix are expressed invariantly, that is, with- 
out reference to coordinates, is essential for their use in 
Sec. ED 



Appendix B: Perturbation theory for degenerate 
modes 

When dealing with the splitting of degenerate eigen- 
values in response to some deformation, as exemplified in 
Fig.[U the question emerges which linear combinations of 
the unperturbed modes must be inserted into Eqs. (|26[) 
and (|2"8)l in order to obtain the different branches. For 
finding these proper linear combinations, we start with 
the following assertion: If u^.\ u± are two nondegen- 
erate solutions to the system (P) - ©, with eigenvalues 
e« ^e<- 2 \ then 

/ d 2 S u [ l ] d n u {2) = f d 2 S u^d n u W = . (Bl) 
Jon Jan 

This is shown by first exploiting Eq. ((3]), and writing 

e {2) d n u {2) - d n u {2) =0 on dn . (B2) 
Multiplying by vP? and integrating, one arrives at 



/ d 2 Su W d n u (2) = 
Jan Jan 



d 2 Su^d n uf 



I d 2 S u { l ] d n u (2) , (B3) 
Jan 
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where Eq. ([2]) has been used. In the same manner one 
also finds 

e« / d 2 S u^d n u W = / d 2 S ufd n u { l ] . (B4) 

JdQ Jdn 

Moreover, one has 

Jan Jn 

= [ d 3 r Vu W ■ VwL 2 ' (B5) 
Jn 

by virtue of Eq. (JTJ) . Hence, we deduce 

/ d 2 S u W d n u {2} = [ d 2 S u [2) d n u { ^ ; (B6) 
Jan Jan 

similarly, 

/ d 2 S u^d n u (2) = f d 2 S ufd n uf . (B7) 
Jan Jan 

Therefore, subtracting Eq. (|B4p from Eq. (|B3|) readily 
yields 

(e( 2 ) - e«) f d 2 S u W d n u {2) = . (B8) 
Jan 

Since e^ 2 ' by assumption, this, together with 

Eq. (|B6j) . was to be demonstrated. 

Next, we consider two such branches of solutions which 
depend on a deformation strength h and are degenerate 
only for vanishing deformation, e' 1 ) = = e for h = 0, 
while ^ e^ 2 ). Since then the assertion (|Blj) holds 
for any h ^ 0, continuity demands that it is also valid 
for h — 0: It is this requirement (|B1I) which singles out 
the proper linear combinations of the degenerate unper- 
turbed modes, when re-tracing the split modes back to 
the point of degeneracy. Focusing now on h = 0, and 
observing 

= f d 3 r (u W Au^ } - (AuL 1} )M- 2) ) . (B9) 
Jn 

virtually the same steps that also lead from Eq. (p~9|) to 
Eq. (|26|) then result in 

[ d 2 Su W d n u^ 
Jan 

= [ d 2 S (u^G™ -e{d n u W )F^) 
Jan 

= {l-e)M{u^\u^) , (BIO) 



where we have introduced the expression 

M(u<V (2) ) = / d 2 S a(V d u W -V a u^+ed n u W d n u^) . 
Jan 

(BH) 

Because the l.h.s. of Eq. (|B10[) vanishes, so does the 
r.h.s. Excluding the particular value e = 1, we deduce 
) = 0. This requirement finally dictates how 
to proceed in the general case: Let us assume that an 
eigenvalue e to the system CU — ([3]) is p-fold degener- 
ate, with eigenmodes v^ k ', where k = l,...,p. Then 
M(i/ fc ) , t>W) is a symmetric px p -matrix and hence pos- 
sesses p orthonormal eigenvectors rf? (where the lower 
index £ refers to the components) with eigenvalues dP\ 
such that 



Y^Miv^M 1 ^) iff = dWrfjp (B12) 
e—i 



and 



E^=««, (B13) 

k=l 



employing the usual Kronecker delta Sij . Therefore, mul- 
tiplying Eq. (|B12p by r/^' and summing over fc, one gets 



J2 M^V^^V^) = dPSij . (B14) 
e,k=i 



Thus, setting 



u O) = Y,4 ] v {k) , (B15) 
fc=i 



we have M{vS x \vP') = for i ^ j, as required. 
This Eq. (|B15[) therefore specifies the desired proper lin- 
ear combinations of the unperturbed, degenerate eigen- 
modes. 
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